arXiv:1502.00889vl [hep-ph] 3 Feb 2015 


Thermodynamics and fluctuations of conserved charges in Hadron Resonance Gas 

model in finite volume 

Abhijit BhattacharyygQ 

Department of Physics, University of Calcutta, 92, A. P. C. Road, Kolkata - 700009, INDIA 

Rajarshi RajQ and Subhasis Samant8@ 

Center for Astroparticle Physics & Space Science, 

Bose Institute, Block-EN, Sector-V, 

Salt Lake, Kolkata-700091, INDIA 
& 

Department of Physics, Bose Institute, 

93/1, A. P. C Road, Kolkata - 700009, INDIA 


Subrata Sui[§ 

Department of Physics, Panihati Mahavidyalaya, 

Barasat Road, Sodepur, Kolkata - 700110, INDIA 

The thermodynamics of hot and dense matter created in heavy-ion collision experiments are 
usually studied as a system of infinite volume. Here we report on possible effects for considering 
a finite system size for such matter in the framework of the Hadron Resonance Gas model. The 
bulk thermodynamic variables as well as the fluctuations of conserved charges are considered. We 
find that the finite size effects are insignificant once the observables are scaled with the respective 
volumes. The only substantial effect is found in the fluctuations of electric charge which may 
therefore be used to extract information about the volume of fireball created in heavy-ion collision 
experiments. 

PACS numbers: 12.38.Mh, 21.60.-n,21.65.Mn,25.75.-q 

Our present day universe contains a significant fraction of matter in hadronic form. In the very early universe — a 
few microseconds after the Big Bang jl] when the temperature was extremely high, the strongly interacting matter is 
expected to have existed in the partonic form. Similar exotic state of matter may exist inside compact stars due to 
extremely high matter density attained by gravitational compression Q ■ For the last few decades various experimental 
efforts are being made to recreate such exotic matter through the collisions of heavy-ions at ultra-relativistic energies. 
Experimental facilities at CERN (France/Switzerland), BNL (USA) and the upcoming facility at GSI (Germany) are 
at the forefront of efforts taken to create these exotic states of matter. One of the major goals in the experiments 
is to study thermodynamic properties of strongly interacting matter at high temperatures and densities. Present 
experimental data as well as lattice QCD simulations seem to indicate a smooth cross over from hadronic to quark 
gluon matter at low density and high temperature SI- At high density and low temperature a first order transition 
is expected @- 113 . 

Usually any thermodynamic study assumes the system volume to be infinite. However the fireball created in the 
relativistic heavy ion collision experiments has a finite spatial volume. The size of such spatial volume critically 
depends on three parameters : the size of the colliding nuclei, the center of mass energy (i /s) and the centrality of 
collisions. Analysis of experimental data could reveal the freeze out volume of the system. One way to carry out such 
an analysis is the study of HBT radii which has been done in Ref. E3- The major finding of this study indicates 
that the freeze out volume increases as the y/s increases and the estimated freeze out volume varies from 2000 fm 3 
to 3000 fm 3 . Another way of estimating the system size is thr oug h the comparison of simulation results with the 
experimental data as done in Ref. [l2j where the UrQMD model Qjl] is used for this purpose. A study of the Pb — Pb 
collisions at different energies and centralities resulted in a freeze-out volume in the range 50 fm 3 to 250 fm 3 . We 
note that these volumes as quoted above, are the freeze out volumes and as one looks back to the early times in the 
evolution of the system, a smaller volume is expected. So both in the quark phase and in the hadron phase (after 
the putative phase transition) finite volume could play an important role in the quantitative measurement of the 
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observables. Thus a thorough study of finite volume effects on the thermodynamic observables in strongly interacting 
matter at high temperature and density is required. Some aspects of such effects have been discussed in the literature 
in the context of QCD inspired models as well as pure gluon theory on the lattice [3®. These studies indicate a 
significant effect of finite volume both in the vacuum and also in the medium. Particularly, in the medium the studies 
in finite volumes have a strong bearing on the location of the critical temperature, critical end point, and related 
thermodynamic variables as discussed by some of us in [35, 36|. Even the existence of critical end point, which in 
turn depends on the order of the phase transition, depends on the volume [35}. In [36] we found that for a two flavour 
system the scaling of ratios of fluctuations with volume is violated around the phase boundary. 

Most of the studies mentioned above are related to the partonic degrees of freedom. However the matter detected 
in the experiments are in the hadronic form and looking at the freeze out volumes one can easily conclude that the 
finite volume effects would be relevant here. So, a finite volume study of the hadronic matter at finite temperature 
and density would be useful. To this end we take recourse to Hadron Resonance Gas (HRG) model [37y| . HRG model 
is based on the Dashen, Ma and Bernstein theorem [13 which shows that a dilute system of strongly interacting 
matter can be described by a gas of free resonances. The attractive part of the hadron interactions is supposed to be 
taken care of by these resonances. 

In this work we plan to study hadronic matter at high temperature and density confined to a finite volume. We 
shall study the thermodynamic properties of hadronic matter and also the fluctuations of various conserved charges. 
We have organised the paper as follows. First the Hadron Resonance Gas model and its finite volume extension is 
introduced. We then discuss the method for calculating fluctuation of different conserved charges. Thereafter we 
present our results and conclude. 

In the existing literature one may find a variety of HRG models. Detailed discussions of the different versions of 
HRG model and some of the recent works using these models may be found in Refs. [37|. 1391454} . In this work we plan 
to discuss the simplest version namely the non-interacting HRG. The basic physics should be independent of which 
version of the model we choose. First we briefly describe the HRG model. 

The grand canonical partition function of a hadron resonance gas [39i . [52} can be written as, 

lnZ ld = ^lnZf, (1) 

i 

where sum is over all the hadrons, id refers to ideal i.e., non-interacting HRG. For the *’th species, 

In Zf = ±^J p 2 dp]n.[l±exp{-(Ei-fiij/T)], (2) 

where Ei = \/p 2 + m 2 is the single particle energy, m* is the mass, V is the volume of the system, <7,; is the degeneracy 
factor and T is the temperature. In the above expression p, t = Bips + Sips + QipQ is the chemical potential and 
Bi, Si, Qi are respectively the baryon number, strangeness and charge of the particle, p's being corresponding chemical 
potentials. The (+) and (—) sign corresponds to fermions and bosons respectively. 

Our purpose is to consider HRG in a system of finite volume. We incorporate this by considering a lower momentum 
cut-off p m i n =7 t/R = A(say) where R is the size of a cubic volume. With this cut-off the partition function for particle 
i becomes 

In Zf = ±^J p 2 dpln[l±exp(-(Ei-m)/T)}, (3) 

From partition function we can calculate various thermodynamic quantities of interest. The partial pressure Pi, the 
energy density Ei and the entropy density .s, can be calculated using the standard definitions, 
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The specific heat at constant volume Cy is given by, 
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Here we consider the grand canonical partition function given in equation J3]) to incorporate the effect of finite 
volume and have incorporated all the hadrons listed in the particle data book [55| up to mass of 3 GeV. 
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FIG. 1: (Color online) Variation of scaled pressure, scaled energy, scaled entropy and scaled specific heat with temperature. 


In Fig. Q]we have shown the scaled pressure i.e., the ratio of pressure calculated for a finite system size R to the 
pressure for infinite volume ( i.e. R = oo) as a function of temperature T. Also shown in the adjacent figures are 
scaled energy, scaled entropy and scaled specific heat. Five different representative system sizes are chosen: R = 2/m, 
R = 4/m, R = 6/m, R = 8/m and R = 10/m. We find a significant volume dependence in all these quantities. At 
extremely low temperature (~ 0.02 GeV) the finite volume effect is strongest. As the temperature increases the scaled 
variables approach unity. This is expected since the lowest lying hadrons are most dominant at low temperatures and 
they feel the finite volume effects the most. With increase in temperature higher mass resonances become important 
and the dependence on finite system size diminishes. 

At around a temperature of 0.1 GeV the scaled quantities for R = 10 fm almost reaches one indicating that the 
results are almost indistinguishable from those at infinite volume. Obviously for smaller volumes the scaled variables 
will reach unity at higher temperatures. We find that by 0.2 GeV of temperature the scaled quantities for system 
sizes down to R = 4.fm reach the value of one. For smaller volumes the system will still be away from the infinite 
volume scenario at T = 0.2 GeV by which the system is supposed to convert to a partonic phase. 

Fluctuations of conserved charges such as net electric charge, baryon number and strangeness have been consid¬ 
ered as probes for hadronization and thermalization of the system created in nuclear collisions [5614581 ] . Moreover, 
fluctuations are expected to show distinctly different behaviour in a hadron resonance gas and a QGP. From the 
grand canonical partition function ( Z ) we take derivatives with respect to chemical potential to get the various 
susceptibilities. The n th order susceptibility is defined as, 
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where p, x is the chemical potential for conserved charge x. We have considered for our purpose x = B (baryon), S 
(strangeness) and Q (electric charge). 

In Fig. [5] we have shown scaled second order and fourth order susceptibilities for conserved charges namely baryon 
number, electric charge and strangeness as function of temperature for pb = 0 for different system sizes R. The 
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FIG. 2: (Color online) Variation of scaled susceptibilities y 2 , X 4 j for baryon (left column), electric charge (middle column) and 
strangeness (right column), with temperature (fi = 0). 


general features of these scaled fluctuations are similar to the previous thermodynamic variables. For all the cases 
we observe a strong volume dependence especially for low temperature. Though the quantitative features of baryon 
number and strangeness are very close to those of the thermodynamic variables discussed earlier, the nature of electric 
charge fluctuations are quite different. We observe that in this case even the systems with volumes larger than 4 fm 
are significantly different from the infinite volume systems. It may be noted that the lightest hadrons - the pions 
contribute to this sector and not to the other charge fluctuations. As expected the ensemble with the lightest particles 
will have the most significant signatures of finite volumes. 

This is further reflected in the observables mean (M), variance (tr), skewness ( S ) and kurtosis ( k ) describing the 
particle distribution. These observables may be combined to relate to the ratios of susceptibilities as, 
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These kind of ratios are important when comparing theoretical predictions with experimental results where the volumes 
may not be estimated reliably and it is expected that the volume factors in the numerator and denominator drop 
out [59], However this assumption is valid when the low momentum scales are less important. The variation of these 
ratios for net proton, net charge and net kaon are shown as a function of the centre of mass collision energy yfs in 
Fig. (0. The parametrization of the temperature and chemical potentials at the freeze-out conditions is taken from 
Ref. [601 - We find that for protons and kaons the ratios are almost independent of volume, but for the net electric 
charge sector a significant volume dependence is present. While the proton and kaon masses are in the range of 0.5 — 1 
GeV, the lowest mass in the electric charge sector is that of the pions with is similar to the lowest momentum scale 
for a finite size system. A similar effect of low momentum scale affecting the results at cross-over has been discussed 
in & 

To summarise, we have studied the thermodynamic properties and fluctuations of hot and dense matter in a 
finite volume. We have used the HRG model for this purpose. The finite volume has a significant effect on the 
thermodynamic properties and also on the fluctuations. But for most observables the effect diminishes once they are 
scaled with the respective volumes. The only significant difference even after scaling with the volume is observed for 
the fluctuations for electric charge. Given that one can relate such volume dependence on the observable in heavy-ion 
collision experiments, it may be possible to extract information about the system volume from the net electric charge 
fluctuations. 
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FIG. 3: (Color online) Energy dependence of a 1 2 3 4 5 6 7 8 9 10 11 /M, Sa and kct 2 for proton (left column), net electric charge (middle column) 
and net kaon (right column). 
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